Noise rectification in quasigeostrophic forced turbulence 
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We study the appearance of large-scale mean motion sustained by stochastic forcing on a rotating fluid (in 
the quasigeostrophic approximation) flowing over topography. As in other noise rectification phenomena, the 
effect requires nonlinearity and absence of detailed balance to occur. By application of an analytical coarse- 
graining procedure we identify the physical mechanism producing such effect: It is a forcing coming from the 
small scales that manifests in a change in the effective viscosity operator and in the effective noise statistical 
properties. 
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Nonlinear interactions can organize random inputs of en- 
ergy into coherent motion. This noise-rectification phe- 
nomenon has been discussed in several contexts ranging from 
biology to physics or engineering [El]. Three ingredients are 
needed to obtain this kind of noise-sustained directed motion: 
nonlinearity, random noise lacking the property of detailed 
balance, and some symmetry-breaking feature establishing a 
preferred direction of motion. 

It has been recently shown numerically [||| that directed 
motion sustained by noise also appears in quasigeostrophic 
two-dimensional fluid flow over topography. The average flow 
at large scales approaches a state highly correlated with topog- 
raphy that disappears if noise or nonlinearity are switched-off. 
The presence of topography provides the symmetry-breaking 
ingredient needed to fix a preferred direction. The small scales 
of the flow follow a more irregular behavior. 

In this Letter we analytically calculate a closed effective 
equation of motion for the large scales of the flow, by coarse- 
graining the small scales. From this effective equation, the 
forcing of the small scales on the large ones (sustained by 
noise and mediated by topography) appears clearly as the re- 
sponsible for the directed currents. The effect is more notice- 
able for increasing nonlinearity, and appears as a renormal- 
ization of the viscosity operator in such a way that it favors 
relaxation towards a state correlated with topography, instead 
of towards rest. This forcing vanishes when noise satisfies 
detailed balance, as in the other kinds of noise-rectification 
phenomena. 

The particular model considered here is the equation de- 
scribing quasigeostrophic forced turbulence. A large amount 
of rotating fluid problems concerning planetary atmospheres 
and oceans involve situations in which vertical velocities are 
small and slaved to the horizontal motion [^P]. Under these 
circumstances flow patterns can be described in terms of two 
horizontal coordinates, the vertical depth of the fluid becom- 
ing a dependent variable. Although the fluid displays many of 
the unique properties of two-dimensional turbulence, some of 
the aspects of three-dimensional dynamics are still essential, 
leading to a quasi two-dimensional dynamics. In particular, 
bottom topography appears explicitly in the equations. This 



kind of quasi-twodimensional dynamics is not only of rele- 
vance to the case of rotating neutral fluids but there is also 
a direct correspondence with drift-wave turbulence in plasma 
physics [|J|]. 

The streamfunction tp(x, t), with x = (x, y), in the quasi- 
geostrophic approximation is governed by the dynamics: M\: 



d\7 2 ip 
dt 



A [V>, V 2 ^ + h] = j/V 4 V' + F , 



(1) 



where v is the viscosity parameter, F(x, t) is any kind of 
relative-vorticity external forcing, and h = fAH/Ho, with 
/ the Coriolis parameter, Hq the mean depth, and AiJ(x) 
the local deviation from the mean depth. A is a bookkeeping 
parameter introduced to allow perturbative expansions in the 
interaction term. The physical case corresponds to A = 1. The 
Poisson bracket or Jacobian is defined as 



_ dAdB _ dB_dA 
dx dy dx dy 



(2) 



Equation ([j]) represents the time evolution of the relative 
vorticity subjected to forcing and dissipation. In the case of 
drift-wave turbulence for a plasma in a strong magnetic field 
applied in the direction perpendicular to x, %jj is related to the 
electrostatic potential, and h — ln(u) c /no), where oj c and 
no are the cyclotron frequency and plasma density respec- 
tively. Eq. (0) is also the limiting case of the more general 
Charney-Hasegawa-Mima equation when the scales are small 
compared to the ion Larmor radius or the barotropic Rossby 
radius [|-|]. 

We now stablish how the dynamics of long-wavelength 
modes in (0), when F is a random forcing, is affected by 
the small scales. Stochastic forcing has been used in fluid 
dynamics problems to model stirring forces [^], wind forc- 
ing [|To|], short scale instabilities thermal noise [ pi] 12| , 
or processes below the resolution of computer models [13Q, 
among others [|l4|]. A useful choice of F, flexible enough to 
model a variety of processes, is to assume F to be Gaussian 
stochastic process with zero mean and correlations given by 

Pk{u))P v (u')) = Dk-y5{k+k')S(LU+Lu'). F k (w) denotes 
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the Fourier transform of F(x, t), k = (k x , k y ), and k = |k|. 
The process is then white in time but has power-law correla- 
tions in space, y = corresponds to white-noise also in space 
and, in the absence of topography, it sustains the Kolmogorov 
spectrum [ p"5||l6| ]. In addition this value of y has been ob- 
served for wind forcing on the Pacific ocean [0. Thermal 
noise corresponds to y = —4 [|l2 16 1. In this case there is a 
fluctuation-dissipation relation between noise and the viscos- 
ity term, so that the fluctuations satisfy detailed balance. 

To obtain the desired large-scale closed equation we have 
applied a coarse-graining procedure to the investigation of the 
dynamics. For our problem it is convenient to use the Fourier 
components of the streamfunction i/j^ or equivalently the rel- 
ative vorticity Ckw = — fc 2- 0kw This variable satisfies: 



^ ^kpq (CpflCqH' + Cpfi/tq) 
p,q,f2,^' 



where the interaction coefficient is: 



.4 



kpq 



{Pxq v -P v qx)p 



k,p+q 



the bare propagator is: 



(3) 



(4) 



(5) 



and the sum is restricted by k = p + q and lu = fi + fi'. 
p = (p x ,p y ), p = |p|, and similar expressions hold for q. 
< k < ko, with fco an upper cut-off. Following the method 
in Ref. jlq] , one can eliminate the modes with k in the 
shell fcoe < fc < fco and substitute their expressions into 
the equations for the remaining low-wavenumber modes £ < 
with < k < k$e~ 6 . To second order in A, the resulting 
equation of motion for the modes is: 



at 



where 



+ A [^<, V 2 V < + h<] = z/V 4 (^< - gh<) + F' , 



X 2 S 2 D(2 + y)5 



32(2tt) 2 ^ 3 



]{\,D,6, v,y) = 



X 2 DS 2 (y + A)S 
16(2tt) 2 ^ 3 



(6) 



(7) 



(8) 



F'(x, t) is an effective noise which turns out to be also a 
Gaussian process with mean value and correlations given by: 



< F'(x,t) >= - 



X 2 DS 2 {A + y)8 



\Q{2ir) 2 v 2 



V 4 h<, 



(9) 



Dk- y S{k + k r )S{uj + uj') (10) 



S 2 is the length of the unit circle: 2ir. Equations (^)-(|l(J) are 
the main result in this Letter. They give the dynamics of long 
wavelength modes tp < . They are valid for small A or, when 
A ~ 1, for small width S of the elimination band. The ef- 
fects of the eliminated short wavelengths on these large scales 
are described in the new structure of the viscosity operator 
and the corrections to the noise term F' . The action of the 
dressed viscosity term V 4 (ip < — gh < ) is no longer to drive 
large scale motion towards rest, but towards a motion state 
<?/i < ) characterized by the existence of flow following the 
isolevels of bottom perturbations h < . This ground state would 
characterize the structure of the mean pattern. The energy in 
this ground state is determined by the function g(A, D, S, v, y) 
which measures the influence of the different terms of the dy- 
namics (nonlinearity, noise, viscosity). Relation (||) shows 
that while nonlinearities jl9| ] and noise increase the energy 
level of the ground state, high values of the viscosity parame- 
ter would imply a reduction of the strength of the ground state 
motion due to damping effect that viscosity exerts over small 
scales. The other mechanism that reinforces the existence of 
average directed motion comes from the fact that the dressed 
noise has got a mean component as a result of the small scale 
elimination. 

A most interesting fact in (^|) and (^) is the presence of the 
factor y + 4. It implies that the tendency to form directed 
currents reverse sign as y crosses the value —4, and that it 
vanishes if y = —4 which is the value for thermal noise satis- 
fying detailed balance. The vanishing of the directed currents, 
obtained here to second order in A, is in fact an exact result 
valid to all orders in the perturbation expansion. This can be 
seen from the exact solution of the Fokker-Planck equation 
associated with Eq.(|l|) for this value of y p0[]. This reflects 
the fact that noise rectification can not occur when detailed 
balance holds. 

As a consistency check we point out that if the term rep- 
resenting the ambient vorticity h is zero, classical results of 
two-dimensional forced turbulence are recovered Jl^]. To 
show that the result implied by the perturbative expressions 
(|^)-(|l0|) is really present for arbitrary A we check the increas- 
ing tendency towards average flow following the topography 
for increasing A: Numerical simulations of (|l]) have been con- 
ducted in a parameter regime of geophysical interest: we take 
/ = lCP 4 ,^ 1 as appropriate for the Coriolis effect at mean 
latitudes on Earth and v = 200m 2 s~ 1 for the viscosity, a 
value usual for the eddy viscosity in ocean models. We use 
the numerical scheme developed in on a grid of 128 x 128 
points with a proper inclusion of the stochastic term [|J] . The 
distance between grid points corresponds to 10 km, so that 
the total system size is L = 1280 km. The amplitude of the 
forcing, D = 1 X 10~ 9 m 2 s~ 3 has been chosen in order to 
obtain final velocities of several centimeters per second. The 
topographic field (shown in Fig. [l]) is randomly generated 
from a specific isotropic power spectrum (Fig. ^) with ran- 
dom phases. The model was run for 6 x 10 5 time steps (corre- 
sponding to 247 years) once a statistically stationary state was 
reached, and some of the results for the mean streamfunction 
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are displayed in Figs. || andQ Currents with a well defined 
average sense appear. Consistently with our analytical results, 
the contour levels of the mean streamfunction follow the to- 
pographic contours more closely the higher the value of A is. 
This is more quantitatively shown in Fig. || where the linear 
correlation coefficient p between the mean field < tp > and 
the underlying topography is plotted as a function of A. A 
spectral analysis of the different resulting fields show that the 
large scales are better adjusted to topography, as well as the 
very small scales where no significant motion is present (Fig. 
|J). Discrepancies are clear for the small but excited scales. 
This can be understood considering that the effect of viscosity 
on these small scales is still to drive the system towards rest. 

For negative values of y noise acts more strongly on the 
small scales, where viscosity damping is more important, so 
that a larger noise intensity is needed to obtain significant 
large-scale directed currents. More important is the reversal 
in the sense of the currents when y < —4. This can be char- 
acterized by the change in sign of p. For example, for y = — 6 
and D = 2 x m 8 s~ 3 , p — —0.6. More detailed results 
will be presented elsewhere [ p0[ ] . 

Concluding, the outcome of this work can be formulated as 
follows: quasigeostrophic flows develop mean patterns in the 
presence of noisy perturbations. As relations (|6|)-(|l0|) show, 
the origin of these patterns is related with nonlinearity and 
lack of detailed balance. Nonlinear terms couple the dynam- 
ics of small scales with the large ones and provide a mech- 
anism to transfer energy from the fluctuating component of 
the spectrum to the mean one. This mean spectral component, 
that is inexistent in purely two-dimensional turbulence [ 22 1, is 
controlled by the shape of the bottom boundary and character- 
izes the structure of the pattern. The existence of these noise- 
sustained structures has a wide range of implications in the 
above mentioned fields of fluid and plasma physics. First be- 
cause it highlights the important and organizing role that noise 
can play in these systems. Secondly, it establishes the need to 
modify not only the value of the parameters (as usually done 
in eddy-viscosity approaches) when performing large eddy 
simulations with insufficient small-scale resolution, but also 
the structure of the equations in a way determined by topogra- 
phy. This last statement has been previously suggested from 
a heuris tic poi nt of view in the context of large-scale ocean 
models [ |23| , |24| . Our results represent a step forward towards 
the justification of such approaches. 
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FIG. 1. Depth contours of a randomly generated bottom topog- 
raphy. Maximum depth is 381.8m and minimum depth —381.8m 
over an average depth of 5000m. Levels are plotted every 63.6m. 
Continuous contours are for positive deviations with respect to the 
mean, whereas dashed contours are for negative ones. 
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FIG. 2. Comparison, as a function of the radial wavenumber 
index Lk 2 /2n, of the power spectra of the bottom topography (solid 
line) and the power spectra of the mean streamfunction obtained for 
A = 0.1 (dotted line), 0.3 (dashed line), 0.6 (dash-dotted line) and 
1 (dash-dot-dot-dot line), y = and D = 10 _9 m 2 s~ 3 . In order 
to carry out the comparison, the fields have been normalized to have 
the same maximum value. 




FIG. 3. Mean streamfunction computed by time averaging 
when a statistically stationary state has been achieved. Continu- 
ous contours denote positive values of the streamfunction, whereas 
dashed contours denote negative ones; A = 0.1, y — 0, and 



D = 10 J mV 



Maximum and minimum values are 1637.7 and 



-1637.7 m /s,and levels are plotted every 272.95 wf/s. 
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FIG. 4. The same as Fig. ^|but for A = 1, maximum and mini- 
mum values are 991.864 and -991.864 m 2 /s, and levels are plotted 
every 165.31 m 2 /s. 
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FIG. 5. Linear correlation coefficient p between the mean 
streamfunction and topographic fields as a function of the interac- 
tion parameter A. y = and D — 10 -9 m 2 s -3 
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